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TITLE OF THE INVENTION 
System and Method For Diagnosing Pathologic Heart Conditions 



CROSS-REFERENCE TO RglsATED APPLICATIONS 
[0001] This application claims the benefkm prior filed co-pending U.S. Provisional 




\ 

Patent Application No. 60/184,3^ filed on February 23, 2000. 



BACKGROUND OF THE INVENTION 
[0002] The present invention relates to a system and method for diagnosing pathologic 
heart conditions based upon heart sound data. 

[0003] Studies have shown that primary care physicians frequently refer patients to 
cardiac specialists on the basis of suspicious heart sounds detected by traditional 
stethoscope auscultation, though a large percentage of these referrals are dismissed by 
cardiologists as having no pathologic condition. The costs, delays, worry, and 
administrative burden resulting from these needless referrals could be reduced if the cues 
that the specialist uses could be incorporated into an algorithm to automatically screen for 
pathologic heart sounds and murmurs. 

[0004] Although attempts have been made to automate screening by auscultation, no 
device is currently available to fulfill this function. Multiple indicators of pathology are 
nonetheless available from heart sounds and can be elicited using certain signal processing 
techniques such as time-frequency analysis. At least one signal of pathology, the systolic 
murmur, can reliably be detected and classified as pathologic using a portable 
electrocardiogram and heart sound measurement unit combined with a time-frequency 
based algorithm. Time-frequency decomposition analysis holds promise for extending 
these results to detection and evaluation of other audible pathologic indicators. 
[0005] In addition, an automatic screening algorithm would be useful for detecting 
pathologic heart murmurs in settings where a trained professional is not always available, 
such as pre sports participation physicals, and examinations performed in remote or 
underserved areas. Furthermore, automated analysis of digitized clinical information such 
as heart sounds could have major implications for health care delivery systems using 
telemedicine. 
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SUMMARY OF THE INVENTION 
[0006] The present invention comprises a time-frequency murmur diagnostic device and 
method. The present invention combines a cardiologist's auscultation expertise, a large 
and growing set of comprehensive heart sound files, and digital signal processing 
algorithms. 

[0007] A method of diagnosing pathologic heart conditions in which a time series of heart 
sounds is filtered and parsed into a sequence of individual heart cycles. A systolic interval 
as well as systolic sub-intervals are identified for each heart cycle. An energy value is 
computed for the systolic sub-interval of one or more heart cycles. The energy value 
computed is proportional to the energy level associated with the filtered series of heart 
sounds. A composite energy value is then computed for the systolic sub-intervals of one 
or more heart cycles and compared to a threshold level in order to distinguish between a 
normal heart and a pathologic heart. 

[0008] The system for diagnosing pathologic heart conditions is comprised of a portable 
computing device that manages data collection and stores data collected from new 
patients, and analyzes data. Also included is a patient data collection unit, communicable 
with the portable computing device, for acquiring and digitizing electro-cardiogram 
(ECG) and heart sound data from a patient. The patient data collection unit is comprised 
of a pair of transducer contact microphones (primary and reference) for obta i ning acoustic 
data. Also included is a pair of acoustic pre-amplifiers connected with the transducers. 
The pre-amplifiers have a passband of 20 Hz to 2 kHz and are used to condition acoustic 
data received from the contact microphones. Variable amplifiers connected with the 
acoustic pre-amplifiers variably amplify the conditioned acoustic data. Moreover, several 
electro-cardiogram (ECG) electrodes connected to an ECG amplifier record ECG data. 
The acoustic and the ECG data are passed to an analog to digital converter connected with 
the variable amplifiers and the ECG amplifier. The data is digitized and sent to the 
computing device for processing by a screening algorithm implemented by the steps 
described in the method above. 

BRIEF DESCRIPTION OF THE DRAWINGS 
[0009] FIGURE 1 illustrates a plot of time-frequency analysis basis functions (wavelets) 
at various scales a. 
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[0010] FIGURILliUustrates a normal heart sound under time-frequency analysis. 
[0011] FIGUR E 3 il lustrates a pathologic systolic heart murmur under time-frequency 
analysis. 

[0012] FIGURE 4 illustrates a block diagram of system hardware. 

[0013] FIGURE 5 illustrates a logic flow diagram of the processes used to diagnose 

pathologic heart conditions. 

[0014] FIGURE 6 illustrates a logic flow diagram for optimizing the parameters of a 
time-frequency screening algorithm. 

DESCRIPTION OF THE PREFERRED EMBODIMENTS 
[0015] Before describing the present invention it is helpful to have a basis for the 
detection of pathological heart conditions. Analysis of a variety of heart sounds and 
corresponding diagnoses from the Johns Hopkins School of Medicine (JHU/SOM) has 
shown that time-frequency analysis is a versatile technique for detecting and classifying 
pathologic heart conditions. Of the available time-frequency techniques, wavelets are a 
useful method for representing heart sound frequency dynamics without creating cross 
term artifacts. Wavelet transforms can be computed for a continuous or discrete set of 
transform variables, depending on the priority for compactness (e.g., data compression 
applications) or ease of visual interpretation (e.g., pattern recognition applications), 
respectively. While compressibility is of interest for future applications, the relative ease 
of interpretation provided by continuous wavelet transforms (CWTs) suggested their use 
with the present invention. While wavelet transform analysis are illustrated herein as a 
preferred method for analyzing heart sound data, Fourier transform analysis may also be 
implemented by one of ordinary skill in the art in order to analyze heart sound data. 
[0016] The CWT of a time series, f(t), is defined as: 



where f and \|/ are both square-integrable, a is a time scaling variable, and b is a time 
translation variable. This can also be written as a convolution: 



\dt (Equation 1) 



W{a 9 b) = f{b)*w a/i (-b) 



where 
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[0017] While it is possible to construct a \\t to yield an optimal peak (compact, high 
amplitude) in W for a given f, it would not be guaranteed to be optimal for a different time 
series g. Since the present invention is to be applied to a variety of heart sound signals, a 
custom wavelet was not implemented. Rather, the alternative was to draw i|/ from a pool 
of wavelets designed to have various advantageous properties. Wavelet transformations 
known as second order "coiflets" were implemented with the present invention. It is 
important to note, however, that other wavelet transformations, including custom 
wavelets, may be implemented without departing from the spirit or scope of the present 
invention. 

[0018] Coiflets of order 2 at various scales, a, are plotted in FIGURE 1. The frequency 
bandpass limits at each scale are given. Due to the equivalence of convolution in the time 
domain and multiplication in the frequency domain equation 1 shows that the Fourier 
transform of \|/ will act as a bandpass filter of the signal f. The bandpass limits in 
FIGURE 1 are the 6 dB passband frequency limits of the Fourier transform of <p . 
Wavelets are constructed so as to maintain a constant ratio of center frequency to 3 dB 
bandwidth (Q), and have a finite duration. Their time-frequency resolution is inherent in 
their design and scale parameters. This is in contrast to Fourier decomposition, which 
uses the infinite time extent sine and cosine functions. Time resolution is not inherent in 
the Fourier transform, but is introduced by the user via windowing the data. Multiple 
Fourier transforms using distinct window intervals would be required to produce the 
constant Q decomposition offered by wavelets. 

[0019] Example cases of a normal heart sound and a pathologic systolic heart murmur 
under wavelet transformation with coiflets are given in FIGURES 2 and 3. These figures 
demonstrate the clear relationship between an audio-based physician's description used in 
auscultation and the visual presentation in time-frequency space. A "harsh pansystolic 
murmur" is a diffuse area of broadscale (i.e., broadband) energy between SI and S2. The 
broadband nature of the sound, which lasts throughout systole, is responsible for 
"harshness." This straightforward representation of a pathologic indicator in time-scale 
space is a promising basis for pattern recognition. 
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[0020] Auscultation of pathologic murmurs is keyed to the following observations, 
according to a study of 222 consecutive patients referred to the Johns Hopkins Pediatric 
Cardiology clinic: 

a. pan-systolic nature of the murmur; 

b. intensity of the murmur > grade 3; 

c. point of maximal murmur intensity at the left upper sternal border (LUSB); 

d. harsh quality of the murmur; 

e. presence of an early or mid-systolic click; or, presence of an abnormal second 
heart sound. 

[0021] The goal was to identify systolic murmurs that are indicative of heart defects, and 
exhibited one or more of the qualities (a)-(d) above. Algorithms may also be implemented 
to detect heart clicks, and split and abnormal S2 sounds for greater diagnostic utility. 
[0022] A system block diagram of the present invention is illustrated in FIGURE 4. The 
system is comprised of two principal elements. One is a patient data collection unit 401. 
The other is a computer processing device 420 including or having access to data storage 
devices. The patient data collection unit performs several functions including obtaining 
heart sound data via a set of contact microphones in the form of transducers. Two 
channels of acoustic data are obtained from a patient using a primary transducer 402 and a 
reference transducer 404. In addition, a set of ECG electrodes 406 are used to obtain 
electro-cardiogram data from the patient. 

[0023] The two contact microphones are each conditioned by a pre-amplifier 408 having a 
passband of 20 Hz to 2 kHz and variable gain amplification stage 410. A set of 
headphones 412 connected to a headphone amplifier 414 can be used to listen to the 
acoustic data gathered from the primary transducer 402 and the reference transducer 404. 
The ECG electrodes 406 feed into an ECG amplifier 416. Outputs from the variable gain 
amplifiers 410 and the ECG amplifier 416 are fed to a analog-to-digital converter where 
the acoustic and ECG signals are digitized and recorded. A 12-bit National Instruments 
PCMCIA analog-to-digital converter is used, for instance, to collect and digitize data at a 
rate (e.g. 8.13 kHz) consistent with the highest data frequencies of interest. 
[0024] Once the signals have been digitized and recorded, the patient data collection unit 
401 forwards the data to a computer processing device 420. The computer processing 
device 420 is typically a laptop computer (due to its compact transportable nature) having 
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adequate data storage capacity. However, the patient data collection unit 401 may be 
connected to other computer processing devices without departing from the spirit or scope 
of the present invention. A computer software program accesses heart sound data that has 
either been collected and forwarded by the patient data collection unit 401, or is resident 
on the laptop computer 420, or can be obtained from another source of heart sound data. 
The computer program applies a screening algorithm to the heart sound data in order to 
determine whether the heart sound data is to be classified as normal or pathologic. Each 
heart sound data file corresponds to a different patient. When the algorithm has operated 
on the heart sound data that has been input, the computer program will display the results 
on a display screen to the doctor, nurse, or technician operating the computer. Results 
indicating a pathologic condition will likely cause the patient to be referred to a 
cardiologist for further examination. Otherwise, a cardiologist referral can be deemed 
unnecessary. 

[0025] FIGURE 5 illustrates the flow of logic and processing that occurs in the various 
elements described in FIGURE 4. Digitized heart sound recordings were collected on 
patients in the Pediatric Cardiology Echocardiography Laboratory in the Johns Hopkins 
Outpatient Center. The recordings were stored in a Heart Sound database 502. Recorded 
heart sounds could then be extracted 504 from the heart sound database 502 and placed 
into an ECG/heart sound data file set 508. Alternatively, ECG and heart sound data could 
be obtained directly from a patient 506 and placed into the ECG/heart sound data file set 
508. The ECG/heart sound data file set 508 serves as the data to be fed to a screening 
algorithm. The purpose of the screening algorithm is to analyze the ECG data and heart 
sound data in order to detect any pathologic anomalies that may be present. Thus the 
system is to be used as a diagnostic aid. In order for the screening algorithm to be applied, 
the data set must first be manipulated. Initially, the ECG data and acoustic heart sound 
data are separated. The ECG data is used to parse the time series into a sequence of 
individual heart cycles via a process that identities the ECG peaks 510. 
[0026] A systolic interval is then identified 512 for each heart cycle. The first and second 
heart sounds are identified either by reference to the heart cycle boundaries or 
acoustically using a passband of 25-140 Hz. In the acoustic method the times of the 
acoustic maxima define systole and diastole. Systole can then be divided into various sub- 
intervals including, but not limited to, the first half of the systolic interval, the total 
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systolic interval, the middle half of the systolic interval, and the last half of the systolic 
interval. A short interval to isolate the first and second heart sounds is factored in. 
Meanwhile, the acoustic heart sound data is passed through a digital bandpass filter 514, in 
this case a second order coiflet CWT transform. Next, a relative energy value (square of 
the wavelet coefficient expressed in dB) for a given wavelet scale and systolic subinterval 
is calculated 516. 

[0027] A composite relative energy value across all included heart cycles is computed. 
The composite energy value can be computed in several manners. One way is to compute 
it as the median of the set of computed energy values for each systolic sub-interval of the 
included heart cycles. A second way is to compute it as the weighted average of the set of 
computed energy values for each systolic sub-interval of the included heart cycles. A third 
way is to compute it as the median relative energy across more than one of the heart cycle 
systolic sub-intervals. A fourth way is to compute it as the weighted average relative 
energy value across more than one of the heart cycle systolic sub-intervals. Those of 
ordinary skill in the art could readily devise other alternative ways in which to compute a 
composite energy value without departing from the spirit or scope of the present invention. 
The specific methods described herein are illustrative and not intended to limit the 
invention to a particular manner for computing a composite energy value. 
[0028] A decision 518 between healthy 520 and pathologic 522 hearts is made on the 
basis of the calculated composite relative energy value being above or below a certain 
threshold. A patient's processed data is saved 524 for optional further technical analysis, 
and added to a database 526. Among the various uses of the database 526 are reviewing 
and improving the algorithm's performance. 

[0029] FIGURE 6 illustrates a logic flow diagram for optimizing the parameters of a 
time-frequency screening algorithm. The diagnostic system and processes described 
above were applied to a set of heart sounds. A set of known pathologic heart sounds and a 
set of known normal heart sounds taken in the first half of each month over a period of 
time were extracted from the heart sound database 602. Files from the latter half of the 
months were preserved as "new" data to test algorithm performance after algorithm 
tuning. Each heart sound data file corresponds to a different patient. The screening 
algorithm is applied 604 to a patient's heart sound data file and the results are recorded 
606. Performance of the screening algorithm over all tested patients for a given threshold 
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was measured by the ratio of called positives to true positives, within a universe of known 
positives (the "sensitivity"), paired with the ratio of called negatives to true negatives, 
within a universe of known negatives (the "specificity"). Sensitivity vs. specificity curves 
for a variety of thresholds and systolic intervals were plotted 607. These curves are called 
"Receiver Operating Characteristic" or "ROC" curves. 

[0030] The test is then repeated 608 on the same heart sound data using a different set of 
wavelet parameters, specifically, the scale (a) and the systolic sub-interval (Ssl) are varied. 
This is done for numerous combinations of scale and systolic sub-interval. Once the heart 
sound data for all patients has been subjected to the screening algorithm and results have 
been recorded for the numerous scale and systolic sub-interval variations, the program 
compares the ROC curves for each wavelet scale/systolic sub-interval combination 610. 
The area beneath each ROC curve is computed and the ROC curve having an area closest 
to "1" is deemed to have the best results. The parameters for that ROC curve are then 
chosen as the optimal parameters to use with the screening algorithm 612. 
[0031] The ROC curves, in the test case, indicated that the most optimal algorithm settings 
used a wavelet scale, a, of 16 and was applied to the middle half systolic sub-interval. The 
best results were typically obtained using a systolic interval centered on systole meaning 
that the midpoint of the systolic interval and the midpoint of the systolic sub-interval are 
the same. Those of ordinary skill in the art could readily adjust the algorithm for different 
wavelet scales and systolic sub-intervals without departing from the spirit or scope of the 
present invention. 

[0032] Using these optimum parameters, the algorithm was then applied to an expanded 
data set that included both half months plus additional data collected while the previously 
described analysis was in progress. Of 143 cases tested, 95 were from normal hearts (with 
and without innocent murmurs) and 48 were from hearts with pathology (i.e., murmur 
grade £ 2). Sensitivity and specificity ratios of 96% were achieved. 
[0033] There are several advantages realized by the present invention. An automatic 
screening algorithm would be useful for detecting pathologic heart murmurs in settings 
where a trained professional is not always available, such as pre sports participation 
physicals, and examinations performed in remote or underserved areas. Furthermore, 
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automated analysis of digitized clinical information such as heart sounds could have major 
implications for health care delivery systems using telemedicine. 
[0034] Automated analysis of heart sound data could also be used by cardiologists to 
quantitatively follow and document changes in severity of certain conditions such as aortic 
stenosis and mitral regurgitation, in which changes in murmur characteristics are known to 
correlate with changes in disease severity. In addition, since data for analysis could 
theoretically be collected using any electronic stethoscope, home care nurses or nurses 
aides could make inexpensive in-home bedside recordings that could be analyzed later to 
monitor changes in certain conditions as reflected in heart sounds. 
[0035] It is to be understood that the present invention illustrated herein is readily 
implementable by those of ordinary skill in the art as a computer program product having 
a medium with computer program(s) embodied thereon. The computer program product 
is capable of being loaded and executed on the appropriate computer processing device(s) 
in order to carry out the method or process steps described. Appropriate computer 
program code in combination with hardware implements many of the elements of the 
present invention. This computer code is typically stored on removable storage media. 
This removable storage media includes, but is not limited to, a diskette, standard CD, 
pocket CD, zip disk, or mini zip disk. Additionally, the computer program code can be 
transferred to the appropriate hardware over some type of data network. 
[0036] The present invention has been described, in part, with reference to flowcharts or 
logic flow diagrams. It will be understood that each block of the flowchart diagrams or 
logic flow diagrams, and combinations of blocks in the flowchart diagrams or logic flow 
diagrams, can be implemented by computer program instructions. 
[0037] These computer program instructions may be loaded onto a general purpose 
computer, special purpose computer, or other programmable data processing apparatus to 
produce a machine, such that the instructions which execute on the computer or other 
programmable data processing apparatus create means for implementing the functions 
specified in the flowchart block or blocks or logic flow diagrams. 
[0038] These computer program instructions may also be stored in a computer-readable 
memory that can direct a computer or other programmable data processing apparatus to 
function in a particular manner, such that the instructions stored in the computer-readable 
memory produce an article of manufacture including instruction means which implement 
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the function specified in the flowchart blocks or logic flow diagrams. The computer 
program instructions may also be loaded onto a computer or other programmable data 
processing apparatus to cause a series of operational steps to be performed on the 
computer or other programmable apparatus to produce a computer implemented process 
such that the instructions which execute on the computer or other programmable apparatus 
provide steps for implementing the functions specified in the flowchart blocks or logic 
flow diagrams. 

[0039] Accordingly, block(s) of flowchart diagrams and/or logic flow diagrams support 
combinations of means for performing the specified functions, combinations of steps for 
performing the specified functions and program instruction means for performing the 
specified functions. It will also be understood that each block of flowchart diagrams 
and/or logic flow diagrams, and combinations of blocks in flowchart diagrams and/or logic 
flow diagrams can be implemented by special purpose hardware-based computer systems 
that perform the specified functions or steps, or combinations of special purpose hardware 
and computer instructions. 

[0040] In the following claims, any means-plus-function clauses are intended to cover the 
structures described herein as performing the recited function and not only structural 
equivalents but also equivalent structures. Therefore, it is to be understood that the 
foregoing is illustrative of the present invention and is not to be construed as limited to the 
specific embodiments disclosed, and that modifications to the disclosed embodiments, as 
well as other embodiments, are intended to be included within the scope of the appended 
claims. The invention is defined by the following claims, with equivalents of the claims 
to be included therein. 
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